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ABSTRACT 

We develop here a new approach for the relativistic modeling of the photons 
moving into a quasi-Minkowskian space-time, where the metric is generated by an 
arbitrary n-body distribution within an isolated Solar System. Our model is built 
on the prescriptions of the theory of General Relativity and leaves the choices of 
the metric, as well as that of the motion of the observer, arbitrary. Adopting a 
quasi-Minkowskian expression of the metric accurate to order (f/c)^, a thorough 
numerical test campaign is conducted to verify correctness and reliability of the 
model equations. The test results show that the model behaves according to 
predictions. Specifically, comparisons to true (simulated) data demonstrate that 
stellar distances are reconstructed up to the specified level of accuracy. Although 
the {v/cY approximation is not always sufficient for its apphcation to future 
astrometric experiments, which require modeling to (f/c)^, this work serves also 
as a natural test-ground for the higher order model, whose formulation is now 
close to completion, and will be presented in a forthcoming paper. 
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1. Introduction 

It is known that the Ught signal carries most of the physical information on the celestial 
objects and the physical field it passed through (see Kopeikin and Schafer (1999); Kopeikin 
and Mashhoon (2002) and references therein). Part of this information, like e.g. the distance 
to the source and its velocity, can be extracted by means of (classical) astrometric techniques, 
but at high-level of accuracy, as demonstrated by ESA space mission Hipparcos, relativistic 
effects cannot be neglected, and this "relativistic astrometry" can be considered part of 
fundamental physics (for example, astrometric measurements are one of the natural choices 
to detect possible deviations from General Relativity generated by a scalar field (f) that 
couples with the metric tensor g to generate gravity (Damour et al. 2002a,b), or they were 
used for a possible attempt at measuring the speed of gravity (Kopeikin 2001; Fomalont and 
Kopeikin 2003). 

The purpose of this paper is to construct a model of the celestial sphere using the 
prescriptions of the General Theory of Relativity in order to take into account the relativistic 
effects suffered by light while propagating through the gravitational field of the Solar System. 
The planned future astrometric experiments set the optimal target accuracy for such models 
to the micro- arcsecond level. To this level, it was shown that the hght propagation will 
be affected not only by the mass of the Sun and of the other planets but also by their 
gravitational quadrupole and their translational and rotational motion (Klioner and Kopeikin 
1992; Kopeikin and Mashhoon 2002). Assuming that the Solar System is isolated and it is 
the only source of gravity, then the accuracy of 0.1 /xarcsec is achieved by considering terms 
of the background metric up to the order of (v/c)^; here c is the velocity of light in vacuum 
and V is an average velocity for energy balance. In the Solar System this is of ~ 10 Km 
sec~^. 

In this sense, the next generation astrometric missions like GAIA, of ESA (Ferryman 
et al. 2001; Ferryman 2003), and SIM, of NASA (Shao 1998; Marr 2003), offer a unique 
opportunity to test our model. On the other hand, the observational error of these satellites 
will be pushed to the microarcsecond level (/i.arcsec) therefore, the need is to implement a 
model of the celestial sphere and of the observables which is accurate to that order. As a ref- 
erence, we recall that Hipparcos (of ESA) , the only previous example of a global astrometric 
mission, had an accuracy of 1 milliarcsecond only. 

The construction of a general relativistic many-body astrometric model with the re- 
quired accuracy is expected to be complicated not only by the mathematical structure of the 
relativistic equations but also by the numerical methods deployed to implement the model 
into software code. Therefore it is crucial to have an efficient strategy for testing the model. 
With this goal in mind, we had already developed a relativistic astrometric model taking as 
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bax^kground geometry the exact (unperturbed) Schwarzschild solution (de Felice et al. 1998, 
2001). That model was used as a basic touchstone for comparison in the construction of 
our many-body model. Moreover, since a model limited to the {v/cf' order of accuracy is 
significantly simpler to handle and test, we decided for a full computer implementation of 
the (f /c)^-model, with the intention to use it as test-ground for the higher order extension, 
which is our ultimate goal. 

Finally, it is worth stressing that our aim was not to add to the theory of light prop- 
agation, which is well known, but to develop a relativistic astrometric model with a new 
theoretical approach that treats the light propagation in a curved space-time in such a way 
that all the possible relativistic perturbations coming from the gravity sources arc naturally 
included in the light-path reconstruction. 

In sections 2 and 3, we define the background geometry in a way compatible with our 
basic physical assumptions. In section 4 we discuss how to handle light propagation through 
the Solar System and provide the differential equations which enable us to reconstruct the 
hght trajectory from a distant star to the observer. Then, in sections 5 and 6, we define the 
observables as the measurements coming from a given satellite and link them to the mathe- 
matical boundary conditions needed to integrate the light trajectory. Finally, in sections 7 
and 8, we present the test campaign used to validate the model. 

In what follows, Greek indices run from to 3 and Latin indices run from 1 to 3. 

2. The space-time metric 

The space-time structure which underlies the development of relativistic astrometry 
must mirror the physical and operational reality experienced by the observer. 

The basic step of this project is to identify the background geometry. Our first assump- 
tion is that the Solar System is isolated; this means that there are no perturbing bodies 
intervening between the emitting stars and the Solar System boundaries. It is clear that this 
assumption may not be fully justified since light rays from distant stars may suffer deflec- 
tions due to microlensing induced by intervening bodies (e.g., as far as the mission GAIA 
is concerned, we expect nearly 1000 microlensing events along the Galactic disk during 5 
years (Dominik and Sahu 2000)); indeed they could generate systematic errors in the data 
reduction of the celestial sphere, however, their number is very small as compared to the 
number of observations (> 10^). Furthermore, we know that most of the stars are in a binary 
system and have proper motion which generates an unsteady gravitational field causing a 
rotation of the celestial reference frame. However, the rotation coefficients have a temporal 
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variation which amounts to 1 //arcsec every 20 years, a period much longer than the hfetime 
of all the future space astrometric missions (Sazhin et al. 1998). So we judge that the above 
cases do not affect the vahdity of the hypothesis that the Solar System is isolated. 

A second assumption is that the Solar System generates a weak gravitational field. We 
can then adopt a quasi-Minkowskian coordinate system so that the space-time metric can 
be written as a perturbation of the Minkowski metric rjap, namely: 

9aP = Vat3 + Kp + 0{h^) (1) 

where the h^p s describe effects generated by the bodies of the Solar System and are small 
in the sense that \haf3\ <^ 1, their spatial variations are of the order of l/ia/sl while their time 
variations are at most of the order {v / c)\has\- Clearly the metric form (1) is preserved under 
gauge transformations of the order of h. The order of magnitude of the correction terms 
entering the hap's, is expressed in terms of powers of {v/c). The hap's are at least of the 
order of {v/cf (Newtonian terms) hence any higher level of accuracy within the model, is 
fixed by the power of {v/c) larger than two that one likes to include in the analysis. 

Having the above considerations in mind, we choose a solution of Einstein's equations 
of the type: 

gap -VaP + Yl + 0{h') (2) 

a 

where now the sum is extended to the bodies of the Solar System. In this approximation, the 
metric tensor (2) has in general a non- vanishing term g^i = 0[{v/c)^] and the non-linearity 
of the gravitational field is confined to terms 0[{v/cy] in /iqo and hij. 

Moreover, in August 2000 the General Assembly of the lAU stated that a solution like 
(2) has to be adopted to define the reference frames and time scales in the Solar System 
(Damour et al. 1991; Brumbcrg and Kopeikin 1989). At the first Post-Newtonian level of 
approximation, the metric tensor nearby any planet of the Solar System takes the form: 

2W 2W^ 



goo = -l + ^-^ + C(c-^) (3) 



9oi = -^ + 0{c-') (4) 
9^J = 5,,(l + ^)+0(c-^) (5) 

with i,j — 1, 2, 3. 

In the equations above, W represents a generalization of the Newtonian potential and Wi 
is a vector potential describing the dynamical contribution to the background geometry by 
the relative motion of the gravitating sources as well as by the peculiarities of their extended 
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structures. The same form (3) to (5) is adopted to describe the metric tensor generated 
by the whole Solar System referring so forth to the Barycentric Celestial Reference System. 
Assuming the space-time to be asymptotically fiat, the potentials W and Wi are a sum of 
integrals containing terms of gravitational mass and mass-current density. Such integrals are 
taken over the support of each body of the Solar System. It should be stressed here that the 
linearized form of the metric, in the weak field and slow motion [v -C c) approximation, is 
the current support for calculating the ephemerides at the Jet Propulsion Laboratory (Soffel 
2000), Cahfornia Institute of Technology. 

A light ray, on its way to the observer from a distant star, first feels the gravitational 
field of the Solar System as a point like mass centered in its barycenter, then, as it gets closer, 
it feels the gravitational perturbations of the individual bodies of the System. To the order of 
{v/cf, the light ray will feel the contributions from the individual mass structures while the 
effects arising from their relative motion and spin would enter terms of higher order. Since a 
time derivative adds a factor of the order of (f /c), time derivatives of the metric coefficients 
generate terms at least of the order of (f/c)^. Our third assumption then is that light rays 
will not feel perturbations of the order 0[(f/c)^] so, within each integration of the light 
trajectories, the Solar System is approximated in our model to a static, non-rotating and 
non-expanding gravitating system. A natural consequence of this assumption is that, along 
with time variations of the metric, we shall also neglect mixed metric coefficients of type /loi- 
Furthermore, in calculating the background metric coefficients which affects the fight signal 
at each spatial point of its trajectory, we need not to consider the spatial location of the 
individual gravitational source at the corresponding retarded time. This does not mean that 
the positions of the perturbing bodies are considered fixed during the whole integration, but 
simply that the ephemeris used do not include the retarded time, and that we neglect, to 
{v/cy, the effects due to their velocities. 



3. Setting the geometry of the celestial sphere 

Let (^',^° = r) be a quasi-Minkowskian coordinate system with respect to which the 
space-time metric takes form (1); as stated, the choice to consider only terms up to {v/cY 
amounts to assume that space-time is static. Hence there exists a time-like Killing vector 
field 7], say, along which the physical properties of the space-time do not change and therefore 
uniquely identifies a time direction. Let the coordinate time r be a parameter along 77 so 
that r]°' = 6q. From the Killing equation f](a;/3) = 0, where semicolon means covariant 
derivative with respect to the given metric and round brackets mean symmetrization, one 
easily deduces that the congruence C^, of Killing lines is, to the required order, vorticity-free. 
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In fact, denoting as P{rj)°'i3 — 5'^ -\- •q'^'qp the tensor operator which projects orthogonally to 
T7, one finds: 

= PinY^^Pir^Yp^np.,, = 0[{v/cf] (6) 

where ujap is the vorticity tensor and square brackets mean anti-symmetrization. A space- 
time which admits a vorticity-free congruence of lines can be foliated, namely it admits a 
family of three-dimensional space- like hypersurfaces S'(^*, r) which are everywhere orthogonal 
to C^. It is always possible to choose a coordinate system such that the surfaces r) 
have equation r =constant; in this case the spatial coordinates can be fixed within each slice 
up to spatial transformations only. We shall term these surfaces S{t). 

Let us now consider the unit vector field u which is everywhere proportional to 77, namely 
■u" = e'^^Q where = (— goo)^^^^ is the normalization factor which assures that UaU" = — 1; 
the associated one-form has components = so it is proportional to the gradient 

of S{t), as expected. Through each point of a slice S{t), for any r, there goes a time-like 
curve, orthogonal to S{t) and having as tangent a vector of the vector field u; the totahty 
of these curves through all the points of S{t) form a non-intersecting family of curves, or a 
congruence Cu, which identifies a physical observer. The property of this observer is to be 
static with respect to the selected coordinate representation, that is the spatial coordinates do 
not vary along its world-lines. The parameter a on the congruence, such that u" = d^" /da, 
is the proper-time of the observer u. We then require that the geometry that each photon 
feels before reaching the target is described by metric (1); moreover we also require that 
within each photon travel time the world-lines of the bodies of the Solar System belong to 
the congruence Cu and in particular the barycenter of the Solar System is fixed at the origin 
of the spatial coordinates on each shce, (Figure 1). With this choice, the observer u will 
be termed locally barycentric. This observer is an essential prerequisite of our relativistic 
astrometric model because at any space-time point and apart from a position-dependent 
rescaling of its time rate, it plays the role of the barycentric observer which is located at the 
origin of the spatial coordinates fixed, as said, at the barycenter of the Solar System. Before 
concluding this section, let us recall here that the constraints on the metric, namely hoi = 
and Crf being Killing and vorticity-free, are gauge invariant only up to the order of {v/cY- 




Fig. 1. — In our model each light ray behaves as if the bodies did not move with respect to 
the center of mass CM. 
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4. The light trajectories 

A photon traveling from a distant star to the observer within the Solar System, would 
see the space-time as a time development of slices of constant r. As stated, we shall treat 
each light trajectory assuming that the bodies of the Solar System were fixed at the spatial 
position they had at the time of observation, say. Evidently, any subsequent light ray will 
be considered updating the positions of the bodies of the Solar System according to their 
actual motion. Let us then prepare the way to a suitable treatment of a light ray. Let 

P{u)a'' = 6j + uauf' (7) 

be the operator which projects orthogonally to u. Because of the unitary condition, the 
parameter a on the trajectories of u is not constant on the slices S{t) but varies differentially 
with the position as a = (7{$,\t). Since the spatial coordinates are constant along the 
unique normal going through the point with those coordinates, the parameter a along it will 
be function of r only; i.e. a = a^i{T). Let us now consider a null geodesic T with tangent 
vector field A;" = d^°'/dX which satisfies the following equations: 

k'^ka = (8) 
— + T-^^k^k" = 0; 

the latter express respectively the null and the geodesic conditions; here A is a real parameter 
on T and F"^ are the connection coefficients of the given metric (Figure 2). Assume that the 
trajectory starts at a point P* on a slice S'(r*) (say) and with spatial coordinates The 
light trajectory will end at the observation place on a slice S{tq) and at a point with spatial 
coordinates ^^q-j. We remember that the origin of the coordinate system is meant to be the 
barycenter of the Solar System. The purpose of this work is to determine ^l, namely the 
coordinates of the star, from a prescribed set of observables. 




Fig. 2. — The point Ci^i''')) belongs to the null geodesic but also to the unique normal to 
the slice S{t) crossing it with a value of the parameter a^i(^\^{T). 
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Since our approximation permits a global foliation of the space-time, we find it more 
convenient for the determination and physical interpretation of the gravitational effects on 
light propagation, to consider the spatial projection of the light ray T on the slice S{to); 
this projection is a curve made of the points of S{tq) which one gets to by moving along 
the unique normal through the point of intersection of the light ray with the slice S{t), for 
any r. This curve, that we denote as T, is smooth and has a tangent vector field i whose 
coordinate components are equal to those of the projection of the tangent k to the null ray 
into the rest frame of the local barycentric observer at each of its points, namely 

r = p{uYf)k^ (9) 

(see Appendix A for a mathematical description of this mapping procedure). The vector I 
physically identifies the local line of sight of the local barycentric observer at each space-time 
point of the hght ray trajectory. Indeed, the knowledge of I is essential to reconstruct the 
whole story of the light ray. The curve T will be naturally parametrized by A. 

Prom u°-Ua. — —1, it follows that: 

r ^k" + u''{u^k'^). (10) 

Clearly it is £^ — showing that lies everywhere on the shce S{to) as expected. Since 
each point of T is the image under the above projection operation of a point of T at time 
T, it is more convenient to label the points of T with the value of the parameter a^i(^x){'^) 
which, as we have already said, uniquely identifies that point on the normal to the slice S{t) 
which contains it. Hence, being 

da^ -{uak'^)dX, (11) 
we define the new tangent vector field: 

df" £°' 



da^^(^x) 



In the same way we denote 



so that 



which implies: 



= r + (14) 
Ma = 1. (15) 
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In what follows we shall denote cr^<(A) as a. We can now write the differential equation which 
is satisfied by the vector field £ Prom (11) and (13), the second of equations (8) writes: 

da da 

+ Vl^ii^ + u^){p + u^)^^. (16) 

In this equation the quantity i^i'^V can be written explicitly in terms of the expan- 
sion ©pCT of Cu (de Felice and Clarke 1990), as: 

Wt^^Ufi = OpJ^r (17) 

where Qp^ — P{u)'^pP{u)^a^(^aU/3) Since the only non vanishing components of the expansion 

are Gjj = {l/2)dohij, the expansion vanishes identically as consequence of the assumption 
to neglect time variations of the metric. From this condition and imposing h^i — equation 
(16) becomes to the required order and after some algebra: 



df"' 1 - 1 - 

- ^?7""9,/ioo = (18) 

If q; = 0, equation (18) leads to ^ = assuring that condition ^° = holds true all 
along the curve T; ii a — k equation (18) gives the set of differential equations that we need 
to integrate to identify the coordinate position of the star: 



- f • 



Here we remind that Latin indexes take values 1, 2, 3. 
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5. Observables and boundary conditions 

Our aim is to determine the coordinate positions of a star, from a prescribed set of 
observables. As first step we shall express the observables in a way which reflects a convenient 
set up for the observer. The latter carries a frame, namely a clock which measures its proper- 
time and a rest-space, which are different from the barycentric proper-time and rest-space, 
respectively given by the parameter (t(t) along the congruence Cu and the space-like slices 
S{t). To determine the boundary conditions needed to integrate equation (19), we need both 
frames; the observer frame allows us to define the measurements and the barycentric one to 
express all the coordinate tensorial components. Each measurement is made at a coordinate 
time To when the observer was at a spatial position with respect to the barycenter given by 
the coordinates ^^g-^. Hence we consider as observables the angles between the direction of 
the incoming photon and the three spatial directions of a frame adapted to the observer. 
These three angles provide the required boundary values for t. Let u' be the vector field 
tangent to the observer's world-fine and let {Aa} (where a — 1,2,3) be a space-fike triad 
carried by the observer. The angles '0(Aa/) ^^^^ incoming light ray forms with each of 
the triad direction, is given by (de Felice and Clarke 1990): 



COSWm - h=Qh = a 20 

where no sum is meant over a and = + '^^'^'a operator which projects into the 

observer's rest-frame (Figure 3). 




Fig. 3. — Fixing the boundary condition at the time of the observation tq, where u' is the 
tangent vector to the world hne of the observer. 
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Prom (14), the above equation can be written more conveniently as 

" U'jil + u'PUp 

where Iq = I{to) and recaUing that 7a^A?Af = 1; this is a matrix equation where the 
unknowns are the photon's spatial directions at the time of observation; they can be 
singled out as 

Ii[u[e'^ - Xia] = ii"Aaa - e'^{u''u^). (22) 

The direction of the light ray, as it is seen from the point of observation, depends on the 
motion of the latter relative to the center of mass. This dependence gives rise to the stellar 
aberration. Did the observer not move with respect to the spatial coordinate grid, namely 
if u'"' — w", then (22) would become IlXia = as expected. Equations (19) and the three 
conditions (22) together with equations (12) and the coordinate positions ^^^^ = C(to) of the 
observer at the time of observation, form a closed system of equations whose solutions are 
the coordinate positions of the star provided one is able to identify the value of the parameter 
a at the emission. The solution of (19), in fact, is of the type C{<^j ^(o)) hence a* = (T(r*) 
which marks the photon emission, is an implicit unknown. The way how to determine cr* 
will be discussed in the following section. 

In order to find the boundary conditions £q for equations (19), we have first to construct 
a tetrad frame adapted to the observer. The tetrad frame {Aq,} forms a system of local 
Cartesian axes, which equips the observer with an instantaneous inertial frame; in fact it is: 

(Ac, I A^) = 77-^. (23) 

Let us relabel for convenience the coordinates as (r, x, y, z); each tetrad vector A^ can 
then be expressed in terms of coordinate components with respect to the coordinate basis, 
as-*-: 



u' = \q = e^'(T,ao + X,5, + Y,ay + ZA) (24) 

Aa = Tado + Xad^ + Yady + Zad, (25) 

where e^' is the normalization function which makes u' unitary and a = 1, 2, 3. Conditions 
(23) must be simultaneously satisfied hence after some algebra we obtain: 

e^' = [- (T^goo + X^gxx + Y,^g^ + Z,^g..)] . (26) 



^The letter s is for "satellite-observer" 
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and 
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(27) 
(28) 
(29) 



where a^h. A general solution of (27-29) is given by: 




T2do + X2d^ + Y2dy + Z2d^ 

ndo + z,d. 



(30) 
(31) 
(32) 
(33) 



where the components are explicitly given in Appendix B. This solution describes the in- 
stantaneous inertial frame of an observer endowed with a general motion in a gravitational 
field described by metric (2). 

It is crucial that our model is consistent, under the same conditions, with the one 
discussed in dc Felice et al. (1998) and de Felice et al. (2001) where the Sun was the only 
source of gravity In order to make comparison easier let us assume that the observer moves 
on a circular orbit around the barycenter of the Sun- Earth system and has, at the time of 
observation, spatial coordinates xq, yo, and zq — 0. The 4- velocity of the observer now reads: 



where cu is the coordinate Keplerian angular velocity of revolution and {xq, y©) are the 
coordinates of the barycenter of the Sun-Earth with respect to the barycenter of the Solar 
System. In this case, the tetrad frame becomes: 




(34) 




(35) 
(36) 
(37) 
(38) 




X,d^ + Y,dy 

T2do + X2d^ + Y2d, 



y 



where: 




(goo (yo - y ) ^gxx (xo - x©) 2gyy ) ] 



1/2 



(39) 




(40) 
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^3 



\f9^x {yo - Vq) 
(yo - yQ)\/-9oo 

(Xq - ^0)^-^00 

1 



(41) 

(42) 
(43) 
(44) 
(45) 



and where we have named: 



We can now solve the system (22) exphcitly, relating the observed quantities to the 
unknowns by means of the observer's comoving frame { }, i.e.: 



1 



1 



3[{xo - Xq)^^^^^ + {yo - yQ)y/g^{ujVE + e^y-ffpo) 
{e'^ioVr, - ^/-goo) 

4(yo ~ ye)V9^y/-^ - (^0 - Xe)V9^{^V^ + e^V-ffoo) 



g^^ (eicvVr, - ^/-goo^ ■ 



(46) 
(47) 
(48) 



6. Identifying the emission time 

Let us assume that the star has no proper motion with respect to the given spatial 
coordinate grid whose origin is at the barycenter of the Solar System. In this case the 
spatial coordinates = ^*(t*) of the star remain fixed with time r and the star's world-line 
is one of the curves of the congruence Cu- If a photon is emitted at r*, say, then its trajectory 
all the way to the observer at (C(o)) ''"0) is mapped into a spatial path on the slice S{to) having 
the boundary conditions fixed as explained in Section 4. Integration along this path leads 
to a solution = ^*((7(t)). Let the same star be observed at a subsequent coordinate time 
'''o = ''"0 + the received photon is emitted at a coordinate time = r* + At and its 
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trajectory is now mapped into a slice S{tq). Evidently, this second observation implies a 
different set of boundary conditions for the integration along the second path hence the 
solution will be a new function = ^*'(cr(r)). Since the spatial coordinates of the star are 
preserved under the mapping, they will be identified by the value of the parameter a{r^) 
such that i'\a{n + At)) = C{a{n)) (see Figure 4). 
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Fig. 4. — Identifying the coordinate position ^* of the star without proper motion and the 
emission time a{r^). 
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A first check of the procedure is shown in Figure 5; here we assume that on the hght 
path acts only the gravitational field of the Sun. Wc considered two stars at coordinate 
distances respectively of 1 and 2 parsecs. Fixing the boundary conditions corresponding to 
observations of the stars in two symmetrically opposite directions with respect to the Sun and 
assuming At equal to six months, we find, in logarithmic units, the points of intersections of 
the two integrated spatial paths. The actual distances are slightly less that 1 and 2 parsecs 
respectively as expected since we calculate proper rather coordinate distances. 
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Fig. 5. — Two examples of a test that deduces the distance of a star from two observations 
of the same objects taken at opposite positions on the orbit of the observer. The first star 
is approximately at 1 parsec, the second at 2 parsecs. Each couple of observations gives two 
trajectories that cross at the expected distance of the objects. 
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7. Testing the model 

The pcrturbative model presented in the last chapter is expected to include all terms to 
the order of (f /c)^. This means that, for a typical velocity in the solar system of ~ 10 km-s""^, 
one can expect deviations from the predictions of exact models of ~ 2 -10""^ arcsec. Equations 
(19) cannot be integrated analytically, so numerical techniques are needed to reconstruct the 
hght path. We then needed a test campaign to check the correctness and accuracy of the 
model. 

We used the implementation of the RADAU integrator given in (Everhart 1985) to 
integrate the system of differential equations for the null geodesic, because it has been 
thoroughly tested, and can be easily implemented to reach very high orders of accuracy. 
Since the geodesic equations include both t = and the system of differential equations 
is of class II according to the definition used in (Everhart 1985), that is 

The tests we devised verify that: 

1. the perturbative model is self- consistent; 

2. the amount of hght defiection caused by each individual body of the Solar System, 
as evaluated in our perturbative model, coincides with that expected for an analytical 
Schwarzschild solution at the same order of accuracy and under the same observational 
conditions; 

3. the model is able to reconstruct stellar distances. 

It is clear from these items that our tests do not consider comparisons with other models 
besides the Schwarzschild solution, like those in Kopeikin and Mashhoon (2002) and Klioner 
(2003); however, we believe that it is still too early for such comparisons. First of all, our 
model is at the {v/c)^ order, so it's not strictly needed to confront it with higher order 
models yet. Moreover, although with the same assumptions some results can be recovered, 
the background mathematical framework is quite different and this is another reason for not 
attempting a detailed comparison. Our model is evolving to the {v/cf order, and this future 
version will be presented in a forthcoming paper and, indeed, compared to the two works 
cited above. 

Before illustrating each specific test, let us make some general remarks on the numerical 
solution of the system of differential equations (19). This system requires a set of six bound- 
ary conditions and, as said in section 5, the natural choice is the three spatial coordinates 
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of the observer at the observation time, and the three components of the vector tangent to 
the spatial hne of sight. 

We are deahng with a numerical problem, therefore we need to define methods to stop 
the integration procedure. 

One possibility is to fix a finite range of integration, the starting point being the position 
of the observer and the ending point that of the star. This approach will be adopted for the 
test on stellar distances. 

The other choice is to stop the integration when the tangent vector to the light trajectory 
becomes constant at the precision level of the computer. This method was used in the light 
deflection test. In particular, when we consider only the gravitational field generated by 
the Sun, the total amount of defiection produced beyond 100 AU is under 0.1 //arcsec level, 
therefore this distance turns out to be a pretty good choice for stopping the integration. 

As a final consideration, we mention that the geodesic equations can be adapted to var- 
ious physical situations; in fact, as explained in section 2, the metric has not been explicited 
the only requirement being that it takes the form: 

gaP^Va/3 + haf3 + 0{h^). (49) 

In our approximation it is 

(r/a/3 + Kf^)rf = 1 + 0{{v/cf), (50) 

so the first step was to verify whether the spatial vectors £ obtained at each step of integration 
satisfied this unitarity condition up to the correct order of accuracy of ~ 10~^^. The results 
show that this always happens along the integration path. 

7.1. Self-consistency test 

The basic test is whether our model satisfies spherical symmetry in the case that the 
Sun is the only source of gravity. 

Therefore, we took the {v/cY approximation of the Schwarzschild metric; the funda- 
mental property of this space-time solution is its spherical symmetry. This means that the 
amount of light defiection measured by an observer in a given position must depend only on 
the angular distance of the celestial object from the gravitational source. 
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Fig. 6. — Geometry taken for the self-consistency test. For each ip there are four stars at 
symmetric positions with respect to the Sun-observer direction. The angle ip goes from ~ 16" 
(i.e. limb-grazing rays) to nearly 180°. The observer is located at the origin of the reference 
frame. 
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We considered an observer at ro = 1 AU and a set of stars placed at different angular 
distances from the Sun. For each ip, we have taken four stars symmetrically positioned 
with respect to the Sun as in Figure 6. 

The results are reported in Table 1 where St/ji are the deflections calculated in the four 
cases. They were calculated to the 0.1 //arcsec level and, as expected, are the same for a 
given angular distance: this is a verification of the physical consistency of the model. 

The way we calculated the hght deflections is discussed below. 
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Table 1. Results for the self-consistency tests, ip is the angular distance from the Sun, 
5-01 . . . 5-04 the deflection obtained in four symmetric cases at the same tfj (only one value is 
reported here for each since they are exactly the same for each 
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7.1.1. Numerical calculation of the light deflection 

The integrator takes as boundary conditions the Cartesian coordinates of the observer 
at the time of observation and the component of the unit vector representing the direction 
of the hne of sight (that is the tangent unit vector to the hght path at the moment of 
observation). The quantities ^ and ^ are recalculated at each step of integration, so the total 
light deflection angle is the angle between the viewing direction at the observation time 
and the direction of the unit vector tangent to the light path at the end of integration 
(Figure 7). 
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Fig. 7. — The deflection angle A is the angle between the initial viewing direction and 
the final tangent unit vector at the end of the integration to the light path C,f but parallely 
propagated along the hght trajectory to the observation point. 
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The most straightforward formula for the computation of an angle between two direction 
would be just 

A = arccos(v|w), (51) 

where v and w are any two given unit tangent vectors. Unfortunately, from a numerical 
point of view, this formula is of little utility, because of the smallness of the deflection angle. 
In fact the numerical accuracy that we need for A is SA ~ 1 /iarcsec, but for angles as small 
as those expected for the relativistic deflections, namely A < 2 arcsec, the accuracy of the 
built-in arccos function drops at the 10^"^ arcsec level. Let us show this with a simple test. 
Consider two points Ri and R2 on the unit celestial sphere with polar coordinates 
and (q;2, S2), and let ai — Si = S2 — and 0:2 variable. Then the angle A between the two 
directions v and w pointing to them is simply given by 0:2. However, application of formula 
(51) leads to the results with poor accuracy shown in Table 2. 
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Table 2. Test for the accuracy of the arccos function for small angles. The angular 
coordinates of two points are given as follows: the first one has ai = di = 0, the second one 
has ^2 = while a2 decreases. The angle A between the two directions v and w pointing 
to Ri and R2 is exactly given by 0:2, but is also calculated taking the arccos of the scalar 
product of the two unit vectors. The results show that when A ~ 10~^ arcsec the accuracy 

of the built-in arccos is less than 1 //arcsec. 
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A solution to this problem is to consider an approximate formula for A which is better 
suitable for our case. From spherical trigonometry the exact formulas for the angular distance 
A and the position angle p of two points on the unit celestial sphere (Figure 8) Ri = (cti, Si) 
and R2 = (Qi2,52) are^ 

sinAsinp = cos52 sin(Q;2 — ai) (52) 
sinAcosp = sin52COs5i— cos52sin5iCOs(Q;2 — CKi). (53) 



^We want to stress here that these formulas require the polar coordinates of a point given by the transfor- 
mations a = arctan(a;/j/), 6 = arcsin(^;/r), and r = \J -\- xp' -\- z^, where x, y, z are the respective Cartesian 
components. 




Fig. 8. — Angular distance A between two points Ri = (ai, 6i) and R2 = (0:2, ^2)- Here P is 
the pole of the celestial sphere and p is the position angle of R2 with respect to Ri. 
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The Taylor series of the sine and cosine functions about x — are 

sinx = x-^ + ^ + 0{x') (54) 

2 4 

cosa; = l- — + ^ + 0(x^). (55) 
2 24 

For X ~ 2 arcsec ~ 10""^ rad, x"^ ~ 10~^° rad ~ 10 /xarcsec and ~ 10~^^ rad ~ 2 • 
10~^ //arcsec we can write V x < 2 arcsec 

sinx = x + C(10~^/xarcsec) (56) 
x^ 

cosx = 1 - — + O(10"Varcsec)). (57) 

Let's now write Q;2 — ai = C 10~^ I'^d ^iid (^2 — (^i = e 10~^ rad. The right-hand 
side of the (53) can be re-written, to the correct accuracy, as 



(2 

~ sin (^2 cos 8x — cos 82 sin (5i ( 1 — — 



sin ^2 cos 5i — cos ^2 sin b\ cos(q;2 — ^i) 

/•2 
2" 

= sin(52 — 5i) -|- -C^ cos ^2 sin b\ 

~ ((^2 — (^1) + ^ COS (^2 sin 5i (58) 

The system of equations (52) and (53) can then be approximated as 

A sinp = (q;2 — Oil) cos(52 (59) 
A cosp = {82 — (^1) + ^ cos ^2 sin ^\ (60) 



so 



{a2 - ai) cos 62 ((^2 - (^1) + 

1 1 

-C^(^2 — (^i) cos ^2 sin 5i + -^'^ cos^ ^2 sin^ 5i 

(as - Oil)' cos^ b2 + (52 - 5i)' + 0(C')- (61) 
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This is the formula that we shall use to compute the angle between any two vectors^. 



•^To the order, we can neglect the difference between cosJi and 00362, in fact 

cos 62 = cos((5i + e) 

~ cos^i ^1 + ^e^^ — esinJi; 

substituting this expression into the Eq.(61) and remembering that (a2 — ai) = C ~ it is easy to show 
that 

(q!2 - aif cos^ 52 = {a2 - aif cos^ + 0{e^). 
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Table 3. Results for the light deflection test. The deflection A^, is computed with a 
formula for the Schwarzschild metric for increasing values of the angular displacement ip of 
the light source from the Sun. This value is compared with that deduced from our model, 
(A). The difference A — A^ goes rapidly under the 0.1 //arcsec level for > 5°, and the 
magnitude of its maximum value is well within that expected for the {v/cY approximation. 
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7.2. The light deflection test 



A formula for the light deflection is given in Misner et al. (1973), for the case of the 
Schwarzschild metric, as 



where Mq is the solar mass in geometrized units, is the distance of the observer from the 
Sun and ip is the angular displacement of a star from the Sun. 

Since (62) is an analytical formula to the same order of (w/c)^, we expect that its 
predictions coincide with those of our model, namely to < 10~^ rad. 

Taking the same stars and computing the hght deflection using the methods seen pre- 
viously, the tests show that the difference between the two predictions is always much less 
than the limit of the approximation. The maximum value (~ 15 //arcsec) is reached for 
limb grazing rays and the difference becomes rapidly less than 0.1 /xarcsec for ip > 5° (see 




(62) 



Table 3). 
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Fig. 9. — The spatial axes of a phase-locked tetrad are such that one is constantly lying in 
the Sun-observer direction, one is always tangent to the circular path of the observer around 
the Sun and the third axes is orthogonal to the orbital plane. The figure shows how these 
axes change when the observer moves from A to B along the fraction of the orbit represented 
with solid line. 
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Table 4. Difference (in //arcsec) between the parallax p* calculated from the exact 
Schwarzschild model and that reconstructed using our {v/cY model (p*). The first column 
reports the exact distance in parsec, the remaining four give (p* — p*) as obtained for each 
of the corresponding s^c) level of accuracy for the impact parameter Tc expressed in 

meters. 
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7.3. Stelleir distances test: comparison with a Schwarzschild model 

Recently we developed a Schwarzschild model for astrometric observations (Vecchiato 
1996; de Felice et al. 1998). In that model we express the components of the vector /c* 
tangent to a null geodesic relative to the spatial axes f, d and of a phase-locked tetrad 
adapted to an observer on a circular orbit around the Sun (Figure 9), namely (de Felice and 
Ussegho-Tomasset 1992) 

er=C0SB(r^fe) = ±- ■ (63) 

■ ■ 3Af0 



ee- = cose(^^,) = ± ^2 ^=== (64) 



3Mq 

To 



e^ = cose(^^^^ 



2 2Mf7 



ro 



A 



To (1 - OOX) 



(65) 



Here Tq is the distance of the observer from the Sun, uj is the coordinate angular velocity 
of the observer, while A and A are two constants of motion of the null geodesic^. The 
constants of motion A and A can be expressed as functions of the impact parameter of the 
hght geodesic with respect to the Sun (rc) and of the angular Schwarzschild coordinates 
of the star (^,0). Finally, Tc can be imphcitly expressed as a function of the complete set 
of Schwarzschild coordinates of the star (r,9,(f)). This means that, given the position of a 
star and of an observer (in Schwarzschild coordinates), we are able to derive the Cartesian 
components of the tangent to the null geodesic at the position of the observer by means 
of an almost completely analytical procedure^. These components allow us to define the 
observables according to Eq.(21). 

In the present model, we can reproduce the same physical situation and the same ob- 
servables by taking hajs as the approximation to the Schwarzschild metric and by adopting 
the tetrad, defined by Eqs.(35)-(38), of an observer on a circular orbit. Hence we can use 
those observables and the position of the observer as the boundary conditions needed to 
integrate backwards the set of differential equations (19). It was this complementarity of the 
two algorithms that was exploited for our test. 



^We have denoted the constants of motion consistently with the notation of the cited work, hence A 
should not be confused with the afRne parameter of the null geodesic mentioned in section 4. 

^The only numerical part of the algorithm is the calculation of rc- 
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We start by giving the "true" position of a star (r*, 9*,(l)*) and two symmetric observers 
(with respect to the Sun) (ri,^i,0i) and (r2,^2,02) in the "Schwarzschild model", from 
which we calculate the tetraedal components of the local light directions for both observers 
using Eqs.(63)-(65). Let's call those components (erj,eej,e^J and (e^.^ , , e<^2 ) • Note that 
(ri, 6i, 0i) and (e^.^, e^^, e<^J represent the boundary conditions needed in our model to inte- 
grate the light trajectory for the first observer, and analogously (r2, O2, 4>2) and (e^j, Gq^-, e^j) 
are the boundary conditions for the second observer. 

The point of intersection of the two paths obtained in this way should be the position 
of the star at (r*, 9* , 0*); however we expect to find slightly different coordinates (r*. 9*^^ 0*). 
Indeed, there are two reasons for this difference: the approximation of our perturbative 
model which is up to (f/c)^, and possible numerical errors. The latter should always be 
maintained below the level of the intrinsic approximation. 

In this test we have taken the two observers at opposite positions on the orbit of the 
Earth around the Sun (i.e. ri = r2 = ro, ^1 = 6*2 = 7r/2 and 02 = 0i + tt). The stars are also 
on the orbital plane {9* — 7r/2) in such a way the two observers are symmetrically placed 
with respect to the Sun-star direction (0* = (0i + 02)/2). Their distances (r*) range from 
1 pc to 10 kpc. 

With this configuration it is easy to calculate the parallax of a star as p* = 1/r* rad 
(the distances being in AU) , and so it is equally easy to determine the numerical accuracy of 
the position in terms of angles by simply subtracting the parallax p* of the exact model from 
the approximate one, p*, reconstructed from the numerical integrations described above. 

The results in Table 4, show that there is an accuracy floor at about 18 //arcsec, as one 
could expect given the intrinsic level of approximation of the model. This table also shows 
that an accuracy on the determination of Tc to ~ 10"^ m at least, is needed to reach the 
theoretical floor for (p* — p*). We stress that in our case, where rc ~ ro — 1.5 • 10^^ m, 
this means that the numerical accuracy needed is £r-c/^c ~ 10~^^, pushing to the edge the 
typical accuracy for a standard double precision number; however to this level the numerical 
contribution to the error is contained to the 0.1 /xarcsec level. 
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Fig. 10. — The geometrical configuration of the deflection test with more than one body is 
such that all the planets are aligned behind the Sun and along the line joining the Sun with 
the observer. 
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Fig. 11. — The upper plot shows the natural logarithm of the partial deflection A-j/^part 
(expressed in arcseconds) of the light ray due to the Sun and the four major planets of 
the Solar System, while the other gives the total deflection. Both are plotted against the 
approximate distance d from the observer in Astronomical Units. All the bodies considered 
are drawn as circles in their approximate locations along the x-axis. 
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8. N-body tests 

The tests described in section 7 consider the Sun as the unique source of gravity. This 
is because in that section we wanted to compare our {v/cY model with an exact one (or at 
least with a model which was already tested at the level of accuracy we require) and the 
obvious choice was the Schwarzschild model. 

However, our model has been built to include an arbitrary number of gravity sources. 
We shall then subject our model with more than one gravitating body to similar tests to 
see whether in this case the results are consistent with expectations. This will also serve as 
the proper testing ground for our forthcoming {v/c)^ model and whatever others would be 
available. Indeed, to the best of our knowledge, we do not know of any extensive numerical 
testing for multi-body relativistic models. 



8.1. Deflection due to a variable number of aligned bodies 

In this test we have calculated the total deflection due to the Sun plus a variable number 
of planets. The geometry is depicted in Figure 10; here the planets are aligned behind the 
Sun with respect to the observer and the light ray was assumed to be grazing the solar limb. 
The results are reported in Table 5 and in Figure 11. 

The flrst row of the table provides the total deflection in four different cases: the flrst 
one gives the result for the Sun plus Jupiter, the second one for the Sun, Jupiter and Saturn, 
and so on. 
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Table 5. Results for the multi-body deflection test. All the numbers are expressed in 
arcseconds. The first row gives the total deflection due to the gravitational pull of the Sun, 
plus that of the planet indicated in the corresponding column heading, and plus that of all 
the other planets on its left. The value of the deflection in the last row, A'^q, is for a hght 
ray grazing the limb of the Sun, i.e. where the distance of closest approach to the center of 

the Sun is ^ Rq = 696 000 km. 
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The second row shows the differences between the deflection due to the Sun (i.e. the 
value Ai/jq) and the amount of total deflection immediately above. This difference gives the 
deflection due to the planets alone. 

The same effect can be approximately estimated calculating the hght deflection due to 
each planet, with the analytical formula for the Schwarzschild solution, and summing them 
up. This is what is represented in the third row. Finally the fourth row simply reports the 
differences among the values in second row and those in third. 

We stress here that the predictions contained in the third row are not well justifled 
from a theoretical point of view. However they can be considered reasonably close to reality; 
hence this and the following tests just tell us that our model behaves "well" with respect to 
a reasonable (but not exact) model under the same physical situation. 
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Fig. 12. — The total deflection due to the Sun and Jupiter was calculated varying the position 
of Jupiter along its orbit. 
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Fig. 13. — When Jupiter is aligned with the Sun and the observer its position on the orbit 
is a = 0°. Its contribution to the total deflection increases approximately to 16.28 mas for 
a ~ 0?31208, where the light path grazes both the Sun and the planet leaving them on the 
same side. For a ~ 0?3226 the light ray still grazes both of them, but passing in the middle, 
so that Jupiter's contribution subtracts to the total deflection. With Jupiter in the range 
0?31208 <a< 0?3226 the photons are stopped by the planet. 
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Table 6. Results for the test of section 8.2. The angular positions of Jupiter are in degrees 

and the deflections in arcseconds. 
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Table 6 — Continued 



Angle 
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-0, 


.0007269 


0.5025188 


1 


.7502596 


-0, 


.0004626 


0.7945559 


1, 


.7505427 


-0, 


.0001795 


0.9401342 


1 


.7505846 


-0, 


.0001376 


1.1236807 


1 


.7506159 


-0, 


.0001063 


1.3762332 


1 


.7506413 


-0, 


.0000809 


1.5891499 


1 


.7506548 


-0, 


.0000674 


1.9463344 


1 


.7506696 


-0, 


.0000526 


2.2474695 


1 


.7506778 


-0, 


.0000444 


2.5127875 


1 


.7506832 


-0, 


.0000390 


3.5539031 


1 


.7506957 


-0, 


.0000265 


7.9518754 


1, 


.7507110 


-0, 


.0000112 


11.2547126 


1 


.7507143 


-0, 


.0000079 


15.9423686 


1 


.7507167 


-0, 


.0000055 


19.5572137 


1 


.7507177 


-0, 


.0000045 


30.0733476 


1 


.7507193 


-0, 


.0000029 


39.7151372 


1 


.7507200 


-0 


,0000022 


54.7655820 


1 


.7507206 


-0, 


.0000016 


67.3801351 


1 


.7507209 


-0, 


.0000013 


78.9125108 


1, 


.7507211 


-0, 


.0000011 


90.0000000 


1 


.7507213 


-0, 


.0000009 
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Table 7. The distance of maximum approach of the hght ray to Jupiter (rmin) is 
calculated in the two extremal cases in which the photons pass nearby the limb of the 
planet. It is always shghtly more than the planet's actual radius R'}^ — 71398 km (Murray 

and Dermott 1999). 



Row 


'"min 




18 


71461 km 


O'.'0 162792 


19 


71474 km 


-0'.'0162771 
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8.2. Deflection of two bodies with variable relative position 

In this test we considered only two bodies, the Sun and Jupiter, and calculated the total 
light deflection experienced by a light ray grazing the solar limb and positioning Jupiter in 
several different places along its orbit. In particular, the planet's position varied in a range 
of 90 degrees, from conjunction to quadrature (Figure 12). 

The results (Table 6 and Figure 13) show that the contribution of Jupiter, as it should 
be, adds to that of the Sun when the two masses are both "on the same side" with respect 
to the light path; their their influences subtract when the light passes between the bodies. 
Moreover, when the light ray grazes the Jovian limb (rows number 18 and 19 of Table 6), the 
contribution to the total deflection due to the planet approaches the predicted theoretical 
value of 16.28 marcsec (de Fehce et al. 2000; Klioner 2000) as shown in Table 7. 

Finally, it is also worth noting that the amount of Jupiter's contribution to the total 
deflection is of the order of 1 //arcsec when the planet is 90 degrees far from the light path, 
again in very good agreement with the theoretical predictions (Khoner 2000). 
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Table 8. 



o+'v Q+%+h e+'^+h+& e+'^+h+&+'S 



distance ~ 1 pc 



d (AU) 


206261.338 


206259.698 


206259.216 


206259.142 


206259.054 


Pi (") 


1.00001775 


1.00002571 


1.00002804 


1.00002840 


1.00002883 


Pi ~ P0 (fJ'dtrcsec) 




-7.95 


-10.29 


-10.65 


-11.07 


^Ipi (jUELTCSGc) 




-8.00 


-10.30 


-10.60 


-11.10 






distance 


~ 10 pc 






d (AU) 


2062283.55 


2062120.14 


2062072.27 


2062064.94 


2062056.3 


Pi (") 


0.10001777 


0.10002569 


0.10002802 


0.10002837 


0.10002879 


Pi — pq (/iarcsec) 




-7.93 


-10.25 


-10.60 


-11.02 


Ai/jj (/iaxcsec) 




-8.00 


-10.30 


-10.60 


-11.00 






distance <■ 


- 100 pc 






d (AU) 


20589921.6 


20573645.4 


20568881.5 


20568152.5 


20567292.3 


Pi (") 


0.01001777 


0.01002569 


0.01002801 


0.01002837 


0.01002879 


Pi ~ P0 (/iarcsec) 




-7.93 


-10.25 


-10.60 


-11.02 


Ai/'i (/uaxcsec) 




-7.99 


-10.30 


-10.60 


-11.00 






distance ~ 


' 1000 pc 






d (AU) 


202664764 


201098829 


200644596 


200575248 


200493473 


Pi (") 


0.00101776 


0.00102569 


0.00102801 


0.00102837 


0.00102879 


Pi — P0 (^iarcsec) 




-7.93 


-10.25 


-10.60 


-11.02 


Ai/)j (/iarcsec) 




-7.99 


-10.30 


-10.60 


-11.00 






distance ~ 


10000 pc 






d (AU) 


1751489760 


1641052370 


1611285280 


1609051510 


1603803810 


Pi (") 


0.00011777 


0.00012569 


0.00012801 


0.00012819 


0.00012861 


Pi — P0 (/iarcsec) 




-7.93 


-10.25 


-10.42 


-10.84 


A'i/'i (/iarcsec) 




-7.99 


-10.30 


-10.60 


-11.00 
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8.3. Stelleir distances revisited: effect of many bodies 

As a final test, we have repeated the experiment of section 7.3, but with an increasing 
number of ahgned planets as in section 8.1 (the sequence of planets was kept the same). 
However, unlike the previous case, in this test we cannot compare our model with any other, 
so this is in fact another test of self-consistency for the (v/c)'^ model. The sense of this 
statement will be clarified below. 
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Fig. 14. — The same initial conditions (i.e. viewing directions) should result in a decreased 
distance of the star when other planets are added. 
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We consider the same boundary conditions as those used for the test with the Schwarzschild 
model. These were unchanged as we added Jupiter, Saturn and the other planets. 

Therefore, the expectation was the the distance to a star decreased as we added more 
planets (and the corresponding parallax increased) because of the increased total deflection 
(Figure 14). 

The results are reported in Table 8. Stellar distances, as before, ranged from 1 pc to 
10 kpc. The first row in each set of 4 is the derived distance, d, in AU (i.e. that obtained 
by the intersection of the two geodesies), the second is the corresponding parallax p — 
206265/d arcsec^. The third row contains the difference between pQ and the parallax of the 
corresponding column. For instance, in the case of Jupiter, the value is pQ — p\ = ^P\- 
Finally, the fourth row provides the contribution A-i/'j, {i = "^,1^,5,^) of the planets to the 
total deflection. This contribution is calculated independently, following the method used in 
section 8.1. 

The self-consistency of the model assures that these last two values coincide up to 
numerical errors. 

The results show that, as we add more planets, the distances decrease confirming our 
qualitative expectations. Moreover, the numerical residuals are Api — Aipi < 10~^ /xarcsec 
(i = %, fi, Z), and this again is compatible with our accuracy for double precision numbers. 

9. Conclusions 

We have developed a general relativistic astrometric n-body model capable of following 
a light ray all the way from the emitting star to an observer (satellite) orbiting around the 
Sun. 

This model has an intrinsic accuracy of ~ 0.1 miUiarcsecond {{v/cY) and its com- 
puter implementation was validated via a thorough test campaign that proved: a) the self- 
consistency of the model; h ) that the amount of light deflection, in the case of a spherical, 
non-rotating, single-body metric, coincides with that produced by an analytical solution; 

c) that, under the same assumptions, our model is able to reconstruct stellar distances; 

d) that in the case of a multi-body configuration, for which we don't have any exact or 



^Since, as in section 8.1, each column gives the results for a given number of planets (i.e., the single Sun 
in first column, the Sun-|- Jupiter in the second, and so on) , we will refer to the quantities in the first row as 
^0, d\, ... and to those in the second row as pq, p\, ... 
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even well-founded analytical approximation, the outcomes for the light deflection and the 
reconstruction of stellar distances are consistent with the results from semi-quantitative 
derivations. 

Finally, this model will be the natural test-bed for the more advanced astrometric model, 
accurate to {v/cf, which will be presented in a forthcoming article. 

Work partially supported by the Itahan Space Agency (ASI) under contracts ASI 
I/R/117/01, and by the Itahan Ministry for Research (MIUR) through the COFIN 2001 
program. 



A. Mathematical description of the mapping procedure 

Here we describe rigorously the mathematical procedure used to obtain the spatial 
projection of the light ray on the shce S{tq). The null geodesic crosses each shce S{t) at a 
point with coordinates = ^*(A(r)); but this point also belongs to the unique normal to the 
slice S{t), crossing it with a value of the parameter a — (t(^'(A), t) = (7^i(;^)(r). (Figure 2). 

Let us now define the one-parameter local diffeomorphism: 

0Aa = </'(V(A)(^o)-a^i(;,)(r)) : T n 5(t) ^ S{t^) (Al) 

which maps each point of the null geodesic T to the point on the slice S{tq) which one gets to 
by moving along the unique normal through the point T{\) fl S{t) , by a parameter distance 
Act = (T^i(;^)(To) — (7^i(A)(T) (Figure 2). Since the spatial coordinates are Lie-transported 
along the normals to the slices, then the points in S{tq), which are images of those on the 
null geodesic under ^ao-, have coordinates (0ao-(^(''^) ^ S'(t)))* = The curve in S{tq) 
which is the image of T under ^ao-, is: 



f = 0A.or; (A2) 

it has tangent vector: 
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Fig. 15. — The curve T will be the image of the null geodesic Tk under the diffeomorphism 

4>Aa- 



-57- 



The map (Al) acts on a 4-dimensional manifold with images in a 3-dimensional one 
hence the coordinates of the target points (Tfl 5'(t)) and those of the image points on S{to) 
are respectively: 

TKro)) = CHt))-t6^. 

Prom this and (A2) it follows that (Figure 15). 

r = 7"^A;^ , (A4) 

hence the curve T is the spatial projection of the null geodesic on the slice S{to) at the time 
of observation. 



B. Explicit expression for the tetrad components 

The following expressions represent the coordinate components of each tetrad vector 
(a= 1, 2, 3) with respect the coordinate basis (a = 0,x,y,z) when the observer moves 
on the orbit which has general barycentric component Xg, Yg, Zg (to be specified): 



Tg 



1 




T- 



= ±- 



'9yy 



' 9xx 

y/gxxX^g + QyyY^ 
Xg a/ QxX 



9yy\/ QxxXg + QyyYg 



\/{9xxXl + gyyY'^)! (goo + g^^Zl) 
^/goo + 9xxX^ + QyyY^ + g-,^Z'^g 
Xg^{gm + gzzZ^)/{gxxX^ + gyyYg^) 

Vgoo + gxxX^ + gyyY^ + g^-^Z'l 
ys\/{goa + gzzZ^^)/{g^^X'^ + gyyY^) 

\/goo + gxxX^g + gyyY^ + g.^Z'l 
Zg^{gxxXl + gyyY^)/ (goo + gzzZ^) 

\/c/oo + gxxXj + gyyY^ + g^.Zj 
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Z3 



= T- 

= 

= 

= ±- 



9qoV 9qo + 9zzZ^ 



900 + gzzZ^ 
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